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HOLOGRAPHIC  PROCESSING  OF  HIGH-SPEED  LIGHTWAVE  SIGNALS  FOR  THE  TIME- 

DIVISION  MULTIPLEXING 


I.  Introduction 


Ultrafast  pulses  contain  a  wide  bandwidth,  with  4  THz  for  a  typical  100  fs  pulse  duration, 
or  10  nm  at  a  center  wavelength  of  850  nm  (~20  nm  at  a  center  wavelength  of  1550  nm).  The 
successful  use  of  this  bandwidth  is  desirable  for  high  throughput  of  coded  data  in  fiber  optic 
communications  channels.  One  method  to  achieve  this  is  to  combine  N  parallel  channels  each  with 
bandwidth  B  into  one  channel  with  bandwidth  NB,  such  as  through  wavelength  division 
multiplexing,  with  the  total  bandwidth  divided  into  N  frequency  intervals  that  are  each  modulated 
independently  over  the  bandwidth  B.  Another  method  is  time  division  multiplexing,  with  the 
modulation  over  the  entire  bandwidth  using  N  time  intervals,  one  for  each  input  channel. 
Unfortunately,  no  modulator  is  available  that  can  operate  at  4  THz.  One  technique  to  shape 
ultrafast  pulses  over  their  full  bandwidth  is  to  use  a  Fourier-domain  pulse  shaper  to  manipulate  the 
pulse  in  the  frequency  domain  to  obtain  the  desired  pulse  characteristics  [1].  This  technique  has 
been  extensively  demonstrated  for  fixed  Fourier  masks  [2]  and  programmable  liquid  crystal 
modulator  Fourier  masks  with  phase  control  [3]  and  amplitude  and  phase  control  [4].  Placement 
of  a  holographic  material  in  the  Fourier  plane  enables  more  complex  and  nonlinear  operations,  such 
as  reversal  of  a  pulse  in  time  [5].  For  dynamic  pulse  shaping,  an  ideal  diffraction  material  has  flat 
amplitude  and  phase  response  over  the  bandwidth  of  the  pulses,  and  a  fast  response  time  for  high 
repetition  rates.  Photorefractive  quantum  wells  are  a  candidate  for  the  dynamic  holographic 
medium  in  a  pulse  shaper  [6,  7].  Diffraction  from  quantum  wells  relies  on  absorption  and  index 
gratings  with  strong  wavelength  dependence.  The  chief  question  investigated  in  this  research  is 
whether  these  strong  dispersion  effects  significantly  distort  the  ultrafast  pulse.  Photorefractive 
quantum  wells  are  evaluated  for  use  in  a  pulse  shaper  by  measuring  the  shape  of  a  diffracted 
femtosecond  pulse. 

Previous  work  in  bulk  photorefractive  crystals  with  two-wave  mixing  of  femtosecond 
pulses  includes  pulse  characterization  [8,  9]  and  pulse  shaping  [10,  11].  Bulk  ferroelectric 
materials  have  the  advantage  of  large  photorefractive  response,  but  the  disadvantage  of  long  grating 
formation  times.  Also,  two-wave  mixing  is  trivially  phase-matched,  but  four-wave  mixing 
requires  tedious  Bragg  matching  of  a  third  beam  in  bulk  materials.  Short  pulses  can  also  suffer 
from  broadening  by  dispersion  when  traversing  a  long  interaction  length  in  a  bulk  crystal.  In 
contrast,  photorefractive  multiple  quantum  wells  have  a  short  grating  formation  time,  and  thin  film 
diffraction  allows  four-wave  mixing  without  a  Bragg  condition.  Quantum  wells  can  also  be 
engineered  in  several  material  systems. 

In  the  nondegenerate  four-wave  mixing  geometry  described  here,  femtosecond  pulses 
diffract  from  a  grating  written  by  an  above-gap  CW  diode  laser.  This  is  in  contrast  to  previous 
work  with  photorefractive  quantum  wells  where  a  near  gap  CW  diode  laser  beam  diffracts  from  a 
grating  written  by  above-gap  femtosecond  pulses  to  find  the  electric-field  correlation  function. 
This  was  accomplished  both  with  direct  interference  [12],  and  interference  in  a  Fourier- transform 
joint  correlator  geometry  (time-to-space  mapping)  [13].  In  both  cases,  photorefractive  quantum 
wells  are  used  to  characterize  the  shape  of  femtosecond  pulses,  but  not  change  their  shape,  by 
obtaining  the  square  of  the  electric-field  correlation  envelope.  The  conceptual  power  of  time-to- 
space  mapping  is  the  correspondence  to  spatial  image  processing  techniques.  These  techniques  are 
exploited  in  spectral  holography. 
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In  the  following  section  II  presents  the  theory  of  diffraction  of  femtosecond  pulses  from 
photorefractive  quantum  wells  and  section  III  describes  measurement  of  pulse  shapes  by  two 
different  interferometric  techniques,  electric-field  cross  correlation  and  spectral  interferometry. 
Section  IV  describes  the  photorefractive  quantum  well  optical  properties  and  the  experimental 
procedures,  and  Section  V,  data  and  analysis. 


II.  Femtosecond  Pulses  and  Photorefractive  Multiple  Quantum  Wells 

In  this  section,  the  temporal  shape  of  a  femtosecond  pulse  will  be  theoretically  calculated 
for  pulses  transmitted  through  and  diffracted  from  a  photorefractive  quantum  well.  The  most 
convenient  method  is  to  Fourier  transform  the  input  pulse  into  the  frequency  domain.  The 
transmission  or  diffraction  of  a  pulse  is  represented  by  a  filter  function,  which  is  multiplied  by  the 
input  electric  field  and  then  inverse  Fourier  transformed  to  obtain  the  temporal  shape  of  the  output 
pulse. 


Femtosecond  pulses  are  emitted  from  a  modelocked  laser  in  a  train  with  repetition  time 
Trep,  and  a  pulse  duration  is  given  as  the  full-width  half-maximum  (FWHM)  tp.  The  center 

frequency  is  C0c>  the  center  wavelength  is  Xc  =  2ncl(dc,  the  propagation  vector  is  k  and  the  unit 
propagation  vector  is  k .  The  electric  field  for  one  pulse  is 


e(x,t)  =  eEofit -k-  x /c)exp  •  x  -  (OJ 


(1) 


where  e  is  the  polarization  vector,  is  the  electric  field  amplitude,  and/(r)  is  the  pulse  shape 
function.  Note  the  sign  convention  in  the  exponential.  In  what  follows,  the  polarization  and 
position  are  known  and  fixed,  so  the  polarization  vector  is  dropped,  and  x  is  set  to  0.  The  electric 
field  simplifies  to 


eit)  =  Eo/(Oexp  (-i©/ 


(2) 


and  the  intensity  is 


=  -lAVG^REp\f^^^\ 

where  Z  is  the  wave  impedance,  and  the  normalization  of  the  pulse  shape  function  is  chosen  to  be 

Jj/(0p^t=l  (4) 


so  that 

—  f  I{t)dt  =  I^YQ  (5) 

^  REPJ-^ 

The  function /(O  is  typically  a  hyperbolic  secant  or  a  gaussian  function  for  modelocked  pulses. 
For  a  gaussian. 
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f{t)  oc  exp 


( 

-  21n2| 
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2'' 

J 


(6) 


where  tp  is  the  FWHM  duration. 

Transformations  between  the  time  and  frequency  domains  are  accomplished  with  the 
Fourier  transforms 


E{ai)  =  F{eit)]  =  £  e(0exp 
e(t)  -  F((o)exp  (-/cot]<ico 


Fourier  ^ 
e{t)  ^  F(to) 


which,  for  gaussian  pulses,  yields  the  relationships 


f  t  f  ^2^ 


7(0  0=  exp 

it 

e(t)  0=  exp  I  -  21n2|-p  j  -  i 


F(co) 


^""P(-41H2 


In2(£^“' 


to. 


iCdct 


Fourier 

<=> 


f 

F(co)  °c  exp  [^-8i5^(?p(C£>  -  03^) 


(7) 


(8) 


The  time  bandwidth  product  for  a  gaussian  is  t^Av  =  tp41n2  /  2TCt^  =  21n2  /  7t  =  0.441  .  If  a 

varying  phase  is  added  to  the  pulse,  the  pulse  may  become  distorted.  Generally,  the  phase  can  be 
expanded  in  a  power  series  of  the  frequency.  The  constant  term  is  trivial.  The  linear  term  is  also 
trivial,  since  it  only  causes  a  shift  in  the  conjugate  variable.  The  lowest  nontrivial  term  is  the 
quadratic  term,  which  causes  chirp.  There  may  also  be  higher  order  terms.  Chirp  may  be  caused 
by  dispersion  of  a  pulse  transmitted  through  a  material.  In  this  case,  the  relationships  are 


7(t)oc 

exp 

it 

e(t)  oc 

exp 


'  41n2ff'^' 


Fourier 


21n2(l  +  /(3)  7-]^- 


I?(co) 


exp 


81n2(^ 


tp(CO-CO^) 

vTTF 


,  2^ 


(9) 
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where  (3  is  a  parameter  characterizing  the  amount  of  chirp.  The  sign  of  |3  is  chosen  so  that  p  >  0 
corresponds  to  positive  chirp  (low  frequencies  ahead  of  high  frequencies),  and  normal  dispersion 
(dn/da  >  0  ).  If  the  chirp  is  caused  by  dispersion,  then  the  bandwidth  remains  constant,  that  is 

tp  /  \/l  +  is  constant,  so  larger  dispersion  results  in  a  longer  pulse  duration.  The  time 

bandwidth  product  for  a  chirped  gaussian  pulse  is  tpAv  =  21n2  \/ 1  +  P^/  n ■  Transform  limited 
pulses  are  the  shortest  possible  pulse  for  a  given  power  spectrum.  Therefore,  increasing 
dispersion  causes  longer  pulses. 

The  transmission  or  diffraction  of  a  pulse  can  be  represented  by  a  complex  filter  function 
Hica)  in  the  frequency  domain.  For  transmission, 

j^co) .  (10) 

Neglecting  Fresnel  reflection  at  the  surfaces  of  the  thin  film  and  the  associated  Fabry-Perot  effects, 
the  filter  fonction  for  transmission  is 

^^(co)  =  exp  ^m(Q))  (11) 

where  n(a))  =  «(«)  +  ia(co)  /  2k  is  the  complex  refractive  index,  k  is  the  propagation  constant  in 
free  space,  and  L  is  the  sample  thickness.  The  transmission  spectrum  is 

r(co)  =  I  p=  exp  (-  a(co)  l]  ,  (12) 

which  can  be  measured. 

In  nondegenerate  four-wave  mixing,  femtosecond  pulses  diffract  from  a  grating  written  in 
the  semiconductor  fhin-film  by  a  CW  laser  with  a  photon  energy  above  the  bandgap  of  the 
semiconductor,  as  shown  in  Fig.  1.  Photorefractive  quantum  wells  operate  by  the 
photogeneration,  transport  and  trapping  of  charge  at  deep  level  defects,  and  an  electro-optic  effect. 
The  interference  of  the  two  writing  beams  creates  a  periodic  intensity  pattern  I(x)  which  induces  a 
refractive  index  pattern 


nix,(i))  =  n(0))  +  M{G))fiKx  +  (\))  (13) 

through  the  photorefractive  effect,  where  A«(co)  =  An(®)  -l-  jAa(co)  /  2k  is  the  complex  change  in 
refractive  index  caused  by  the  electro-optic  effect,  which  implicitly  depends  on  the  applied  electric 
field  as  well  as  beam  ratio,  fringe  spacing,  and  material  parameters  involving  transport.  The 
grating  vector  is  K,  and  ([)  is  the  phase  offset  between  the  intensity  pattern  and  the  index  pattern. 
The  function /(iTx  +  (j))  is  a  periodic  function  with  period  2%,  and  is  generally  not  a  sinusoid  due  to 
the  nonlinearities  in  the  transport  and  the  electro-optic  effect.  The  writing  beams  completely 
determine  the  refractive  index  in  Eq.  13,  where  the  frequency  co  is  for  the  probe  beam.  The 
spectral  dependence  in  Eq  13  is  determined  by  the  transmission  and  electro-optic  spectra  of  the 
thin-film  stmcture.  The  photon  energy  of  the  writing  beams  is  implicitly  included  in  K. 

In  the  thin  grating  limit,  the  spatially  varying  phase  shift  diffracts  the  incident  pulse  into 
many  orders  without  requiring  phase-matching.  The  strongest  diffraction  is  from  the  lowest 
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harmonic  (fundamental)  grating,  so  f(Kx  +  (]))  is  taken  to  be  cos{Kx  +  (j))  [14],  The  total  electric- 
field  at  the  exit  face  of  the  material  is 

Eall(.^)  =  £’;^(co)exp  ^m(o))  kljexp  |An((o)  kLcos  +  (t))j  .  (14) 

Using  the  relation 

exp  ^izcos  (8)j  =  /""/^(zjexp  (/w5)  (15) 

and  the  previous  definition  for 

OO  /  N 

kljexp  Um[Kx  +  (j)] j  (16) 

which  contains  all  diffracted  orders.  The  first-order  diffracted  electric-field  is  for  m  =  1, 

Eiyicd)  =  E j((£))  /7j|An((n)  ^fjexp  -i-  (|)jj  .  (17) 


Using  the  approximation  J J^z)  =  /  2^m\  for  small  z,  and  neglecting  the  phases,  which  are 

unimportant,  the  electric  field  is  written  in  terms  of  a  filter  function  £’£>((0)  =  £';^(co)//£)(co) ,  where 


/7£)(co)  =  //£(®)i^Art(co)  klj 

=  7/7-(co)ifAn(co)  kL  +  /Aa(co)  ^ 


and  the  power  spectrum  for  the  diffracted  pulse  is 


(18) 


££,(03)  =  Ei^iOi)  7/o(®) 


2  2 1  T 

E[}^{&)  Hj{&)  ^  A«(o3)  kL -f /Aa(03)  ^ 


£,^(co)|V(co)i 


(a«(co)  kLj^  -t-  Aa(03) 


(19) 


and  the  input  diffraction  efficiency  is  r\  ^^(co)  =  |  ££,(co)  ^  / 1  £/yv(o))  | 

diffraction  efficiency  is  rj  out(^)  =  Ej){(i))  |^  / 1  Ejid))  |^  =  |  ^  / 1  Elj-ioi) 


In  the  time  domain,  the  electric  field  is  ^£,(0  =  £  ^{££,(03)}  .  The  shape  of  the  pulse  is 

/£,(0  =  I££)(co)I  /  2Z  .  Unfortunately,  it  is  difficult  to  directly  measure  the  temporal  profile  of  an 
ultrashort  pulse.  Therefore,  several  indirect  techniques  must  be  used,  including  linear  electric-field 
interference  of  the  signal  pulse  with  a  reference  pulse,  and  nonlinear  correlation  of  the  intensity 
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profile  with  a  reference  pulse.  With  a  known  reference  pulse,  the  signal  pulse  can  be  fully 
characterized. 


Spectra  and 
Spectral  Interference 


Fig.  1.  Experimental  geometry  for  diffracting  pulses  from  the  photorefractive  quantum  wells.  The 
shape  of  the  diffracted  pulses  is  inferred  from  electric-field  cross-correlation  as  well  as  spectral  and 
spectral  interferometry  methods.  An  above-gap  CW  diode  laser  writes  the  grating  for 
nondegenerate  four-wave  mixing.  The  CW  beam  is  shown  as  the  dashed  line. 


Three  linear  techniques  can  be  used  to  determine  the  shape  of  the  diffracted  pulse,  as 
shown  in  Fig.  1.  These  techniques  include  finding  the  electric-field  cross-correlation  of  the 
diffracted  pulses  with  the  reference  pulse,  finding  the  spectra  of  the  diffracted  and  reference  pulses, 
and  finding  the  spectral  interference  of  the  diffracted  pulse  and  the  reference  pulse.  The 
interference  of  the  signal  pulse  (the  transmitted  pulse  or  the  diffracted  pulse)  and  the  reference 
pulse  (the  input  pulse)  result  in  a  total  electric  field 


+  Ejft  X) 

=  Esfs{t)Q^V  (-  i^st)  +  ERfR^t  -  x)exp  (-  /®^(t  -  x)) 


(20) 


where  x  is  the  delay  of  the  reference  pulse  relative  to  the  signal  pulse.  The  reference  pulse  is  the 
input  pulse,  and  the  intensity  of  the  interference  is 


=  ^1  Esit)  p  +  ^1  ER{t  -  X)  1^  +lRe  [EsmRit  - 


=  I  sit)  +  iRit  -  X)  +  fsit)fRit  -  x)exp  (-  iAcotjexp  (-  im^x)  j 
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where  A®  -  ®s  -  cOr.  In  the  following,  the  frequency  difference  is 
A®  =  0.  A  standard  silicon  photodetector  cannot  respond  faster  than  10 
than  100  fs.  Therefore,  the  detected  intensity  is  averaged 


incorporated  into  fs,  and 
ns,  which  is  much  longer 


hoTAd'^)  -  J  ^ TOTAd^^'^)dt 

=  ^5  +  4  +  2\/7^Re  exp(-/®^T)r  fs(t)fl(t-T)dt 

J—OO  I 


=  / 


TOTAU 


1  +, 


mRe  exp(-/®^T)Y^^(x) 


(22) 


=  I 


TOTAU 


1  +  m 


V 


cos  arg(Y5,/f(T))  -  ®^i; 


where  m  -  2\/ 1  si r  /  {Is  +  Ir)  is  the  modulation  depth  of  the  interference,  and 


ys,R(d)  =  £  fdi)f*R{t  -  T.)dt 


(23) 


is  the  correlation  function,  which  contains  all  the  useful  information  about  the  interference 
gaussian  signal  and  a  gaussian  reference  pulse,  when  both  are  transform  limited, 


For  a 


75, /?(^) 

ocexp  [-41n2-T 

1 

(24) 


which  has  a  width  +  ,  and  reduces  to  x^^^  =  for  an  autocorrelation.  The 

correlation  function  can  also  be  found  for  the  case  where  the  signal  pulse  is  a  chirped  gaussian  but 
the  analytic  expression  is  complicated,  and  therefore  not  directly  useful.  Generally,  the  width  of 
the  correlation  will  be  smaller  than  expected  for  a  transform  limited  pulse  of  the  same  duration  due 
to  the  reduced  coherence  between  the  signal  and  reference  pulses.  In  the  special  case  where  signal 
IS  a  chirped  gaussian  with  the  same  power  spectrum  as  the  reference,  the  duration  of  the  signal 

pulse  is  ts  =  Wl  and  the  correlation  width  is  Xs^r  =  +  (|3/2)2  .  This  is  the  case 

where  the  pulse  is  transmitted  through  a  transparent  material  with  dispersion. 

A  connection  between  the  time  domain  and  the  frequency  domain  can  be  made  through  the 
Fourier  transform  of  the  correlation  function  ^ 


1s,r{'^) 


Fourier 


F5(w)F^(®) 


and  in  the  case  of  autocorrelation  reduces  to 


(25) 


1r,r{^)  |Fr(®)Poc|£^(co) 


(26) 
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which  is  just  the  power  spectrum.  This  is  useful  for  checking  consistency  between  time  domain 
electric-field  autocorrelation  and  the  power  spectrum. 

In  the  frequency  domain,  there  is  interference  of  the  spectra  [15,  16].  Following  Ref. 
[16],  the  electric-field  in  the  frequency  domain  is 

£rorAL(®>'^)  =  (icot)  (27) 

=  E 5jF’5(©  -  (O/j)  +  -  ©/j)exp  (icox] 

and  the  detected  intensity  from  spectral  inteferometry  is 

2 

f  =  1  r  +  ERi(d)exp  (itox) 

+  2Re  |£5(o))Ei(co)exp 

=  I  Esioi)  1^  + 1  Ef^m  f  +  2|  Es{(£>)  II  £/?(co)  |cos  -  ©x) 

WntesTre  t  ^“quantify  is  the  sum  of  the  individual 

interference  spectrum 

S(©,x)  =  I  Etotal(^’^^  f  - 1  ^5(®)  f  - 1  Eni^)  | 

=  2Re  [Es{(ii)E  R(©)exp  (-  i©x)  j 
=  2  E 5  Re  f  F^C©  -  ©/^)F«(©  -  ©i?)exp  (-  i©x)] 


which  has  a  functional 

an  envelope  with  magnitude  and  the  phase  of  the 

as  the  electric  field  correlation). 


(30) 


x(i:)  =  2Re  |^x(x)  exp  j 

=  x(x)  exp  (-/co^x)  +  x*(x)  exp  (/co^x) 
f|x(x)|  =  X((o  -  (oj  +  X*(-  CO  -  (oj 
e(co)F{x(x)}=X((o-o)J 
F~  ^|0(co)F|x(x)||  =  Jc(x)  exp  [-/co^x] 

where  0((o)  is  the  step  function  that  selects  the  positive  frequency  values.  For  a  function  Z(co) 
(such  as  the  spectral  interference), 


X(co)  =  2Re  |x(co)  exp  [ztox)  j 

=  X((o)  exp  (/cox)  +  X*(co)  exp  [-/cox) 

F-  i{X(co)}  =  x(t  -x)  +  x\-t-x)  (31) 

0(OF-i{X(co)}  =  Jc(/-x) 
f|0(OF-  '  {X(co)}|  =  X(co)  exp  (/cox) 

In  both  cases,  the  transform  of  the  real  function  is  the  sum  of  the  desired  function  shifted  to  the 
right  and  its  complex  conjugate  shifted  to  the  left.  The  overall  operation  is  valid  when  the  function 
and  its  conjugate  do  not  overlap  at  the  origin,  so  that  the  step  function  selects  all  of  the  desired 
function  and  none  of  its  conjugate.  For  the  electric-field  cross-correlation  using  Eq.  30,  the 
condition  is  that  the  bandwidth  of  the  pulse  is  much  smaller  than  the  center  frequency.  This  is 
equivalent  to  a  slowly  varying  envelope  approximation  where  there  are  many  oscillations  under  the 
envelope  of  the  function.  This  is  true  for  100  fs  pulses  in  the  infrared,  but  breaks  down  for  10  fs 
pulses.  Similarly,  for  the  spectral  interference  using  Eq.  31,  the  condition  is  that  the  delay  must  be 
much  larger  than  the  pulse  width,  again  so  that  there  are  many  oscillations  under  the  envelope.  In 
this  case,  the  delay  can  be  chosen  to  satisfy  this  condition. 

Both  interferometric  techniques  depend  on  the  characteristics  of  the  reference  pulse. 
Spectral  interference  of  two  pulses  provides  the  relative  phase  of  the  pulses.  If  the  reference  pulse 
is  close  to  transform  limited,  the  phase  difference  is  the  phase  of  the  signal  pulse.  The  combination 
of  second-harmonic  autocorrelation  and  measurement  of  the  spectmm  of  the  reference  pulse  can  be 
used  to  estimate  how  close  the  pulse  is  to  transform  limited.  When  the  reference  pulse  is  not 
transform  limited,  a  technique  such  as  frequency  resolved  optical  gating  (EROG)  [17]  must  be 
used  to  fully  characterize  the  reference  pulse.  Standard  techniques  are  sufficient  for  this  work 
because  the  reference  pulses  are  assumed  to  be  well  behaved  and  the  frequency  difference  between 
the  signal  pulse  and  the  reference  pulse  is  the  quantity  of  interest. 

Second  harmonic  generation  correlation  depends  on  the  intensity  of  the  pulse.  For 
background  free  second  harmonic  correlation,  the  second  harmonic  output  intensity  is 
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(32) 


/&)  =  (22o>('.t)) 

=  (/s(()/„(t-T)) 

where  <>  denotes  time  average  over  a  repetition  period.  The  interferometric  second  harmonic 
generation  correlation  depends  on  the  relative  phase  of  the  two  pulses.  If  the  amplitudes  are  taken 
to  be  equal,  the  second  harmonic  intensity  (using  the  total  electric  field  from  Eq.  20)  is 


oc(|E203aT)|') 

OC  ^1  j  ^ 

^1  1"^  + 1  fkit  -  "t)  I  +  4|  f^t)  I  I  fj^{t  -  x)  I 

+  .  (33) 

+  -  x)exp  (-  i©,x) 

+  2fs(t)Mt  -  X)  f*/{t  -  x)exp  (-  /®,x)\ 


The  measured  intensity  at  the  peak  is  eight  times  the  value  for  large  pulse  separations.  When  the 
interferometric  intensity  is  averaged  over  one  interference  fringe, 

Jo  \  /  (34) 

=  /|  +  4  +  4/f^(x) 

which  is  the  background  free  result  with  a  background  due  to  the  second  harmonic  of  each  pulse. 
For  gaussian  pulses. 


/f^(x)ocexp  -41n2-2^ 


(35) 


which  has  a  width  of  xf^  =  sjtl  +  tl  ,  and  reduces  to  x®  J  =  /2  for  autocorrelation. 


TV.  Experiment 

The  multiple  quantum  well  structure  used  in  our  experiment  was  grown  by  molecular  beam 
epitaxy  in  a  Varian  Gen  II  chamber.  The  growth  began  with  0.5  p,m  GaAs  on  a  semi-insulating 
GaAs  substrate,  followed  by  500  nm  of  Alo.5Gao.5As,  20  nm  of  AlAs,  and  10  nm  of  GaAs  was 
grown  for  substrate  removal  etches.  After  250  nm  of  Alo.1Gao.9As,  85  periods  of  7.5  nm  GaAs 
wells  and  10  nm  of  Alo  iGao  9 As  barriers  were  deposited,  followed  by  150  nm  of  Alo.1Gao.9As. 
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After  growth,  the  samples  were  proton  implanted  with  a  dose  of  10^2  at  160  keV  and  5  x 
10^^  cm'2  at  80  keV  to  create  deep  level  defects  and  render  the  material  semi-insulating  [18,  19] 
and  photorefractive  [6,  20].  The  implanted  wafer  was  cleaved  into  3  mm  x  3  mm  pieces.  To 
remove  the  substrate  for  transmission  experiments,  the  substrate  was  lapped  to  100  |im,  and  then 
etched  with  a  1:19  ammonium  hydroxide/hydrogen  peroxide  etch  to  the  50%  A1  stop  etch  layer 
[21].  Then  an  HF  acid  etch  removed  the  50%  A1  layer  and  the  AlAs  layer  [22].  Gold  contacts  are 
evaporated  directly  onto  the  low  bandgap  GaAs  layer  with  a  1  mm  aperture  across  which  a  field  is 
applied  in  the  plane  of  the  quantum  wells,  perpendicular  to  the  growth  direction. 

The  transmission  and  electro-optic  effect  of  the  photorefractive  quantum  wells  are  shown  in 
Fig.  2.  The  data  are  measured  using  a  1000  W  halogen  source  filtered  by  a  0.85  m  double 
spectrometer.  There  is  a  peak  absorption  due  to  the  quantum-confined  heavy-hole  exciton  of 
approximately  10"^  cm'i  for  the  1.9  |im  thickness  of  the  quantum  well  region.  The  change  in 
absorption  and  the  change  in  the  refractive  index  spectra  due  to  the  electro-optic  effect  are  shown  in 
Fig.  2(b).  The  electro-optic  effect  with  a  field  applied  in  the  plane  of  the  quantum  wells  is  due  to 
the  field  ionization  of  the  excitons,  which  broadens  both  the  light-hole  exciton  absorption  feature 
and  the  heavy-hole  absorption  feature,  resulting  in  a  characteristic  double  dip.  The  differential 
transmission  is  measured  by  applying  a  sinusoid  that  varies  between  0  V  and  400  V  with  a 
frequency  of  277  Hz.  The  change  in  absorption  is  calculated  using 

Aa(F,^)  =  -  i-ln  +  1  j  (36) 

where  L  is  the  thickness  and  F  is  the  applied  field.  The  change  in  refractive  index  was  calculated 
using  the  Kramers-Kronig  transformation 


An(co)  =  fpf  -^^ds  (37) 

with  a  maximum  absorption  change  of  1000  cm'*  and  a  maximum  refractive  index  change  of 
0.005.  The  CW  input  diffraction  efficiency  spectrum  shown  in  Fig.  3(a)  is  calculated  directly  from 
the  transmission  and  electro-optic  spectra  using  Eq.  19.  Since  the  heavy-hole  exciton  oscillator 
strength  is  much  larger  than  light-hole  exciton  oscillator  strength,  the  peak  is  almost  entirely  due  to 
the  heavy  hole  exciton,  with  a  maximum  input  diffraction  efficiency  of  lO'^.  The  spectrum  of  the 
diffracted  phase  shown  in  Fig.  3(b)  is  calculated  only  from  the  changes  in  absorption  and  refractive 
index,  and  not  from  the  background  index  spectrum,  which  is  not  known  in  detail.  Despite  the 
oscillation  in  the  electro-optic  spectrum,  the  phase  is  relatively  flat.  That  is,  the  change  in  refractive 
index  spectrum  is  approximately  the  same  as  for  the  change  in  absorption,  but  shifted  by  a  constant 
phase,  because  of  the  relationship  through  the  Kramers-Kronig  transform. 
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Fig.  2.  Optical  properties  of  the  photorefractive  quantum  wells:  (a)  transmission  spectrum;  (b) 
spectra  for  change  in  absorption  and  change  in  refractive  index  due  to  electro-optic  effect,  with  4 
kV/cm  applied  field.  The  refractive  index  change  is  calculated  from  the  change  in  absorption  using 
the  Kramers-Kronig  transform. 
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Fig.  3.  Calculated  diffraction  spectra  of  the  photorefractive  quantum  wells:  (a)  input  diffraction 
efficiency;  (b)  and  diffraction  phase.  The  diffraction  efficiency  is  calculated  using  both  the 
transmission  spectrum  and  electro-optic  effect.  The  phase  is  only  calculated  from  the  electro-optic 
effect.  Note  that  the  phase  is  flat  near  850  nm. 
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The  femtosecond  pulse  is  diffracted  from  the  photorefractive  quantum  well  using  non¬ 
degenerate  four-wave  mixing,  as  shown  in  Fig.  1.  The  femtosecond  pulse  source  is  a  Clark  NJA- 
4  modelocked  Ti-Sapphire  laser  (with  100  MHz  repetition  rate),  pumped  by  a  Coherent  Innova- 
310  Argon  Ion  laser.  A  685  nm  diode  laser  is  the  source  for  the  writing  beams,  with  a  photon 
energy  above  the  bandgap  of  the  quantum  well  barriers.  The  fringe  spacing  of  the  interference 
pattern  is  10  |im,  which  is  far  above  the  resolution  limit  of  3  pim  and  below  the  limit  where 
scattered  light  from  the  transmitted  beam  overlaps  with  the  diffraction.  The  incident  intensity  from 
the  writing  beams  is  2  mW  total  with  a  beam  ratio  of  unity,  while  there  is  2  mW  in  the  probe 
pulsed  beam,  which  results  in  a  total  intensity  of  approximately  40  mW/cm^.  To  obtain  the  largest 
diffraction  efficiency,  the  writing  beams  should  be  much  stronger  than  the  probe,  but  diffracted 
power  is  maximized  when  the  ratio  of  the  writing  beam  power  to  the  probe  power  is  near  unity. 
An  advantage  of  photorefractive  quantum  wells  is  their  ultra-low  sensitivity,  where  the  gratings  are 
fully  developed  for  a  saturation  intensity  of  10  [xW/cm^.  There  is  a  trade-off  between  response 
time  and  intensity,  with  shorter  response  times  for  higher  intensities.  At  an  intensity  of 
40  mW/cm2,  the  response  time  is  about  10  ^is.  Over  the  duration  of  one  100  fs  pulse,  the 
photorefractive  grating  remains  essentially  static.  To  obtain  diffraction,  a  field  is  applied  by 
applying  400V  DC  over  a  1  mm  gap  for  an  average  field  of  4  kV/cm,  and  the  writing  beams  are 

chopped  at  1.0  kHz. 

The  shape  of  the  pulse  is  measured  using  the  linear  techniques  described  above,  which  use 
an  undistorted  reference  pulse.  The  delay  line  for  the  reference  pulse  uses  a  hollow  retroreflector 
mounted  on  a  stepper-motor  driven  translation  stage  with  a  0.1  |xm  step  size.  A  biased  p-i-n 
silicon  photodetector  is  used  for  electric-field  correlation.  Spectra  are  measured  with  a  0.275  m 
spectrometer  with  a  CCD  for  capture.  Nonlinear  correlation  is  accomplished  with  0.5  mm  thick 
BBO  with  a  photomultiplier  tube  to  detect  the  second  harmonic.  Both  field  and  second-harmonic 
correlation  scans  are  performed  by  chopping  the  writing  beams  and  measuring  the  signal  with  a 
lockin  amplifier.  For  correlation  scans,  the  envelope  function  can  be  used  to  estimate  the  signal 
pulse  shape,  given  the  properties  of  the  reference  pulse.  However,  the  combination  of  the 
uncertainty  of  the  translation  stage  and  drift  of  the  relative  path  lengths  as  the  translation  stage  is 
scanned  makes  phase  information  unreliable  for  correlation  scans.  There  is  no  drift  between 
adjacent  points  in  the  spectra  because  the  CCD  captures  all  points  simultaneously. 


V.  Data  and  Analysis 

The  characteristics  of  the  reference  pulse  are  important  for  interferometric  and  correlation 
techniques.  Pulses  from  the  laser  used  in  this  experiment  are  typically  100  fs  in  duration  and  are 
nearly  transform  limited  gaussian  pulses.  The  reference  pulse  electric  field  autocorrelation  is 
shown  in  Fig.  4(a),  and  the  interferometric  second  harmonic  autocorrelation  is  shown  in  big.  4(b). 
The  sign  convention  is  that  positive  delays  correspond  to  the  reference  pulse  arriving  after  the 
sisnal  pulse  at  the  detector.  Therefore,  the  leading  edge  of  the  pulse  is  on  the  left,  so  the  pulse 
appears  to  be  moving  from  right  to  left.  With  a  100  nm  step  size  for  the  translation  stage  and 
approximately  850  nm  center  frequency,  there  are  over  four  points  per  interference  fringe.  The 
heavy  line  in  Fig.  4(a)  is  the  envelope  of  the  correlation  function  jR^Rix)  calculated  using  Eq.  30, 
which  has  a  width  of  213  fs.  The  power  spectrum  is  shown  in  Fig.  5(a)  with  a  width  of  8  8  nm, 
and  the  average  of  the  second  harmonic  autocorrelation  is  shown  m  Fig.  5(b)  with  a  width  108  ts. 
Both  curves  in  Fig.  5  are  fitted  with  gaussian  functions  to  show  that  the  pulse  is  approximately  a 
gaussian  shape.  Electric  field  autocorrelation  of  the  reference  and  the  Fourier  transform  of  the 
spectrum  match  well,  as  shown  in  Fig.  6.  These  curves  should  match  according  to  Eq.  26,  which 
verifies  that  the  spectral  data  is  consistent  with  the  correlation  data.  The  time-bandwidth  product  is 
calculated  directly  from  the  spectrum 
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Fig.  4.  Reference  pulse  auto-correlation  vs.  delay:  (a)  electric  field  auto-correlation  of  reference 
pulse;  (b)  interferometric  second  harmonic  autocorrelation  of  reference  pulse.  In  (a),  the  heavy  line 
is  the  calculated  envelope. 
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Fig.  5.  Gaussian  fits  to  reference  pulse;  (a)  gaussian  fit  to  reference  pulse  spectrum;  (b)  gaussian 
fit  to  reference  pulse  second  harmonic  autocorrelation,  which  is  the  average  of  the  data  in  Fig. 
4(b).  Both  curves  cire  well  fit  by  3.  g3ussi3n.  The  reference  pulse  is  nearly  transform  limited. 
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Fig.  6.  Electric  field  autocorrelation  of  reference  pulse  and  inverse  Fourier  transform  of  reference 
pulse  power  spectrum.  The  agreement  between  the  two  curves  indicates  consistency  between  the 
correlation  and  spectral  measurements. 

using  tpV  =  tpCAX  /  ,  where  c  =  300  nm/fs,  and  tp  is  calculated  from  the  second  harmonic 

autocorrelation  width  using  Eq.  30.  The  time-bandwidth  product  for  the  reference  pulse  is  0.41, 
which  is  consistent  with  a  transform  limited  gaussian  pulse. 

The  absorption  and  refractive  index  spectra  distort  the  transmitted  pulse.  The  electric  field 
cross-correlation  of  the  transmitted  and  reference  pulses  is  shown  in  Fig.  7(a),  and  the 
interferometric  cross-correlation  is  shown  in  Fig.  7(b).  The  transmitted  pulse  has  a  tail  on  the 
following  edge  of  the  pulse.  The  power  spectrum  in  is  shown  in  Fig.  8(a)  and  the  phase  spectrum 
is  shown  in  Fig.8(b),  calculated  from  the  spectral  interference  of  the  transmitted  pulse  with  the 
reference  pulse.  The  tail  is  due  to  either  the  amplitude  or  phase  of  the  transmission  function 

In  Fig.  9,  the  power  spectrum  data  is  compared  with  the  power  spectrum  calculated  from 
the  product  of  the  transmission  in  Fig.  2(a)  and  from  the  power  spectrum  of  the  reference  pulse. 
In  Fig.  10,  the  electric-field  cross-correlation  data  is  compared  with  caleulations  from  spectral 
measurements.  The  calculation  using  both  the  amplitude  and  the  phase  nearly  matches  the  data, 
while  the  ealculation  without  the  phase  fails  to  produce  a  tail.  Therefore,  photorefractive  quantum 
wells  distort  transmitted  pulses  due  to  both  absorption  and  refractive  index  variation  over  the 
bandwidth  of  the  pulse. 


17 


B 

o 

3 

c 

o 


■iS 

"S 

u 

b 

O 

U 


e 

o 

B 

u 

a 

a 


■a 


(b) 


Delay  (fs) 


9605 14n  thru  nlxc 


Fig.  7.  Cross-correlation  of  transmitted  pulse  with  reference  pulse  vs.  reference  pulse  delay:  (a) 
electric  field  cross-correlation;  (b)  interferometric  second  harmonic  cross-correlation.  The 
transmitted  pulse  shape  has  gained  a  tail  on  the  following  edge  of  the  pulse. 
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Fig.  8.  Spectrum  of  transmitted  pulse:  (a)  transmitted  and  reference  pulse  spectra;  (b)  transmitted 
pulse  phase  measured  by  spectral  interference.  There  is  a  notch  in  the  transmitted  pulse  spectrum 
due  to  the  heavy  hole  exciton  absorption  feature. 
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Fig.  9.  Spectrum  of  transmitted  pulse  and  calculated  spectrum. 


Fig.  10.  Electric  field  cross-correlation  of  transmitted  pulse.  The  first  calculation  uses  the 
amplitude  and  phase  from  the  spectral  measurements,  while  the  second  calculation  uses  only  the 
amplitude,  with  the  phase  set  to  zero.  The  phase  is  necessary  to  describe  the  tail  on  the  correlation. 
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The  peak  of  the  diffraction  efficiency  spectrum  in  Fig.  3(a)  broadens  the  diffracted  pulse  by 
reducing  its  bandwidth.  The  electric-field  cross-correlation  of  the  diffracted  and  reference  pulses  is 
shown  in  Fig.  11.  The  calculated  envelope  has  a  width  of  433  fs,  more  than  twice  as  long  as  the 
reference  field  autocorrelation  width  of  213  fs.  Assuming  the  diffracted  pulse  is  a  transform 
limited  gaussian,  Eq.  24  can  be  used  to  estimate  that  the  diffracted  pulse  duration  is 


which  is  285  fs,  where  the  reference  pulse  duration  is  the  same  as  used  above.  The  increased 
duration  of  the  diffracted  pulse  compared  to  the  reference  pulse  is  in  general  due  to  a  combination 
of  the  decreased  bandwidth  of  the  pulse,  and  the  variation  of  the  refractive  index  over  the 
bandwidth,  that  is,  chirping  (or  higher  order)  changes  in  the  phase.  In  Fig.  12(a),  the  spectra  of 
the  reference  and  diffracted  pulse  are  shown.  The  diffraction  spectrum  appears  to  also  be  nearly 
gaussian,  and  the  spectrum  and  the  estimate  of  the  pulse  length  are  used  to  calculate  a  time- 
bandwidth  product  of  0.43,  also  very  close  the  that  for  a  gaussian.  The  phase  in  Fig.  12(b)  is 
nearly  flat,  confirming  that  the  pulse  is  nearly  transform  limited. 

The  spectrum  of  the  diffracted  pulse  is  calculated  from  the  electro-optic  data  and  compared 
with  the  data  in  Fig.  13.  The  difference  between  the  data  and  the  calculation  is  likely  due  to  the 
difference  between  the  modulation  depth  of  the  space-charge  field  during  photorefractive  mixing 
and  the  modulation  of  the  field  used  for  differential  transmission.  Even  with  a  beam  ratio  of  unity 
in  the  writing  beams,  the  effective  modulation  of  the  space-charge  field  will  be  limited  by  the 
transport.  If  the  effective  space-charge  modulation  were  known,  then  a  precise  comparison  could 
be  made.  Some  error  may  also  be  due  to  nonuniform  distribution  of  the  field  between  the  contacts 
[23] ,  but  this  effect  should  be  small  since  the  incident  intensity  is  high. 

The  electric-field  correlation  of  the  diffracted  and  reference  pulses  is  calculated  from  the 
spectral  data  using  Eq.  18  and  compared  with  the  data  in  Fig.  14.  The  calculation  with  and  without 
the  phase  are  nearly  identical,  indicating  that  the  phase  variation  of  the  diffracted  pulse  is  small,  so 
the  pulse  is  nearly  transform  limited. 

In  conclusion,  pulses  diffracted  from  photorefractive  quantum  wells  are  broadened,  but 
remain  nearly  transform  limited.  The  broadening  is  due  to  the  bandwidth  of  the  diffraction 
efficiency  spectrum,  which  is  determined  by  the  electro-optic  properties  of  the  quantum  wells. 
Since  the  phase  of  the  diffracted  pulses  remain  essentially  flat,  therefore,  the  chief  limitation  to  the 
diffracted  pulsewidth  is  the  photorefractive  bandwidth,  not  the  dispersion  of  the  index  gratings. 
This  indicates  that  photorefractive  quantum  well  devices  will  be  useful  for  pulse  shaping.  Future 
work  includes  the  design  of  structures  with  wider  bandwidth  through  bandgap  engineering  of  the 
quantum  well  structure  [24] .  To  obtain  more  diffracted  power,  asymmetric  Fabry-Perot  stmctures 
[25],  and  perpendicular  geometry  p-i-n  structures  [26]. are  being  explored.  Finally,  photorefractive 
materials  that  operate  at  a  wavelength  of  1.5  |Lim  [27]  are  required  for  pulse  shaping  compatible 
with  fiber-optic  communications  technology. 
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Fig.  12.  Spectra  of  the  diffracted  pulse;  (a)  power  spectrum  of  the  diffracted  pulse  compared  to 
reference  pulse;  (b)  phase  of  the  diffracted  pulse.  The  diffracted  pulse  is  also  nearly  a  gaussian 
shape. 
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Fig.  13.  Spectrum  of  the  diffracted  compared  with  calculation  from  transmission  and  electro-optic 
data:  (a)  linear  scale;  (b)  log  scale. 
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Fig.  14.  Electric-field  cross  correlation  of  diffracted  pulse  with  reference  pulse.  The  calculation  is 
from  spectral  data.  The  small  difference  between  the  calculation  with  and  without  phase  indicates 
the  pulse  is  nearly  transform  limited. 
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